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We show that the large sequential decay corrections obtained by Ono et al nucl-ex/0507018 , is 
in contradiction with both the other dynamical and statistical model calculations carried out for 
the same systems and energy. On the other hand, the conclusion of Shetty et al. [Phys. Rev. C 
70, 011601R (2004)], that the experimental data favors Gogny— AS interaction (obtained assuming 
a significantly smaller sequential decay effects), is consistent with several other independent studies. 

PACS numbers: 25.70.Pq, 25.70.Mn, 26.50.+X 



In Ref. U, we compared the experimentally deter- 
mined isoscaling parameter a, to those of the dynamical 
AMD calculation of Ono et al. Q . The comparison was 
carried out without the sequential decay corrections. It 
was shown that the experimental data agrees better with 
the choice of Gogny-AS interaction. The conclusion was 
based on the assumption that the sequential decay ef- 
fect is small and less than 15% for the systems studied 
(in fact, a correction of 15% results in even better agree- 
ment with Gogny-AS interaction). The same data was 
also compared with a statistical model calculation in Ref. 
[3, and the sequential decay effect was observed to be 
small. The magnitude of this effect was found to be sim- 
ilar to those assumed in the dynamical model comparison 
in Ref. Q. 

Ono et al., in their comment show that the effect of 
sequential decay in the dynamical AMD model is signif- 
icantly large (about 50 %), and therefore the conclusion 
of Shetty et al., 0, that fragment data show preference 
for the Gogny-AS interaction, cannot be justified. The 
conclusions of this comment are based on the premise 
that the secondary de-excitation effect for the dynam- 
ical model is large and significantly different from the 
statistical model calculations. Recently, another dynam- 
ical calculation using IQMD model, has shown that no 
significant difference exists between the primary and the 
secondary a 0, thereby indicating that the effect of se- 
quential decay is very small. The observed small differ- 
ence between the primary and the secondary a is in direct 
contradiction with the large sequential decay corrections 
reported by Ono et al., 4]. The sequential decay effect 
from the IQMD calculation was carried out for the same 
systems and beam energy as studied by Ono et al. Q 
using the AMD model. 

The stark contrast in the sequential decay effect be- 
tween two different dynamical calculations, as well the 
statistical calculation, for the same systems and energy, 
is not clear at this moment. If one assumes that the 



50 



-J 40 



cn 

0) 
c 

UJ 

0) 

in 



30 



20 - 



10 - 



Heiselberg et ol, (Ref. 9) 
Danielewicz (Ref. 8) 
Chen etat (Ref. 14) 
Piekarewicz etat (Ref. 12) 
Shetty et al, (Ref, 




0.5 



1.5 



P/Po 



*On leave from Institute for Nuclear Research, 117312 Moscow, 
Russia. 



FIG. 1: Comparison of the density dependence of the symme- 
try energy obtained from various different studies. See text 
for details. 



sequential decay is determined by the excitation energy, 
charge (Z) and mass (A) of the fragments, and not by 
the process that leads to these fragments, the calcula- 
tions of Ono et al. are also in conflict with their own 
sequential decay calculation from statistical models 
It appears that the isoscaling parameter, like in all sta- 
tistical calculations, remains a robust observable against 
sequential decay even in dynamical calculations. While 
this needs to be verified further, we show in the following 
section that by assuming a negligibly small sequential de- 
cay (as established in various statistical calculations and 
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TABLE I: Parametrized form of the density dependence of the symmetry energy obtained from various independent studies. 

Reference Paramctrization Studies 

Heiselberg et al. [9] 32.0(p/po)°'^'' Variational calculation 

Danielewicz et al. [8] 31(33) (p/po)° '^^'° '^''' BE, skin, isospin analog states 

Chen et al. 31.6(p/po)^"^ Isospin difussion 

Piekarewicz et al. \12} 32.7(p/po)° '^'' Giant resonances 

Shetty et al. [1] 31.6(p/po)'^'^^ Isotopic distribution 



recently shown in dynamical IQMD calculation), the re- 
sults of Shetty et al. , that the data agrees better with 
the choice of Gogny-AS interaction, are consistent with 
several other independent studies. 

Figure 1 shows the density dependence of the sym- 
metry energy obtained from the work of Shetty et al. 
along with those obtained from various other stud- 
ies. The shaded region in the figure correspond to the 
result of Danielewicz obtained by constraining the 
binding energy, neutron skin thickness and isospin ana- 
log state in finite nuclei. The dotted curve correspond 
to the paramctrization adopted by Heiselberg et al. joj 
in their studies on neutron stars, and obtained by fitting 
the predictions of the variational calculations of Akmal 
et al. 0]. The dashed curve correspond to those ob- 
tained recently from an accurately calibrated relativistic 
mean field interaction, for describing the Giant Monopole 
Resonance (GMR) in ^°Zt and ^osph and the IVGDR in 
208pb, by Piekarewicz a/. |Tl|l3|. The solid curve 
correspond to those obtained from Gogny-AS interac- 
tion that explain the results of Shetty et al. (Ij, as- 
suming a negligibly small sequential decay effect. The 
dot-dashed curve correspond to those used for explaining 
the isospin difussion data of NSCL-MSU [3 by Chen et 
al. |lj|- This dependence has now been modified to 
include isospin dependence of the in-medium nucleon- 
nuclcon cross-section, and is in close agreement with that 



of Shetty et al. . The parameterized form of the den- 
sity dependence of the symmetry energy obtained from 
all these studies are as shown in table I. 

The close agreement between various studies shown in 
figure 1, suggests that the assumptions made and the 
conclusion arrived at in Shetty et al. are probably 
correct. 

Finally, Ono et al. Q, argues that the method of " 
extrapolating " the asymmetry of the primary fragments 
as done in Shetty et al. I], is problematic. They base 
their argument on estimation of / ~ 0.75 for ^"^Zn -|- 

Au, ^^Zn + 92Mo and ^^Zn + ^^Ni reactions. These 
systems are highly asymmetric and have charge and mass 
nearly twice that of the systems studied by Ono et al^, 
and Shetty et al. 0, 0| . The systems studied in Shetty 
et al. dllg {^^Fe + ^^Fe, ^^Fe + ^^Ni, ^^Ni + ^^Ni, ^OAr 
-I- sspe, 40 Ar + ^^Ni and 4°Ca + ^^Ni), and Ono et al. 
(4"Ca -1- 40Ca, ''^Ca -I- ^^Ca and ^oCa + *^°Ca), are nearly 
similar, and it is therefore natural to interpolate as done 
in fig. 3 of Shetty et al. 0. Ono et al.'s assertion that the 
calculations are strictly for Ca -I- Ca reactions is contrary 
to the spirit of any theoretical model calculations. 
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